Microscopic origin of isotropic non-Heisenberg behavior in highly correlated systems 
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We have reanalyzed the microscopic origin of the isotropic deviations that are observed from 
the energy spacings predicted by the HDVV Hamiltonian. Usually, a biquadratic spin operator is 
added to the HDVV Hamiltonian to account for such deviations. It is shown here that this operator 
cannot describe the effect of the excited atomic non-Hund states which brought the most important 
contribution to the deviations. For systems containing more than two magnetic centers, non-Hund 
states cause additional interactions that are of the same order of magnitude as the biquadratic 
exchange and should have significant effects on the macroscopic properties of extended systems. 

PACS numbers: 71.10.-w,75.10.-b,71.15-m,71.27+a 
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The magnetic interactions between S=l/2 sites can ac- 
curately be parametrized with the standard Heisenbcrg- 
Dirac-van Vleck (HDVV) Hamiltonian fl. Usually this 
Hamiltonian is extrapolated to systems with higher spin 
moments and the interaction between such magnetic sites 
gives rise to an energy spectrum with a regular spac- 
ing between the different levels, the so-called Lande pat- 
tern: E(S — 1) — E(S) = SJ, where J parametrizes the 
strength of the magnetic coupling between magnetic cen- 
ters. However, in some cases significant deviations from 
this regular pattern are observed and extra terms must 
be added to the HDVV Hamiltonian. One of the most 
commonly applied extensions of the HDVV Hamiltonian 
is the addition of the biquadratic exchange term: 



(1) 



where (ij) are couples of interacting sites, J e ^ is the 
effective bilinear exchange and X e " the biquadratic ex- 
change. Numerous theoretical and experimental studies 
have established that the biquadratic interaction signifi- 
cantly affects the magnetic properties of both ferromag- 
nets and antiferromagnets [2, l3|. For instance, the fer- 
romagnetic phase transition [J] changes character from 
first-order to second-order for a critical value . The 
spontaneous magnetization, the exchange energy and the 
spin-correlation function exhibit discontinuous jumps at 
\ e ff = and unstable behavior for > . One 
may also quote that from spin wave theory, both ferro- 
magnetic and antiferromagnetic spin structures change 
abruptly to canted ones [5(] for a critical value of X e ^ . 
More recently, the phase diagram of the S=l model given 
in equation[T]has been precisely studied in triangular lat- 
tice in a magnetic field, with emphasis on the quadrupo- 
lar phases Q as well as in spin one chains where the open 
question of the existence of a ferroquadrupolar phase 
between the dimerized and the ferromagnetic phases is 
adressedQ- 

The theoretical explanation of the appearance of a bi- 



quadratic exchange was initially given by Anderson [8| 
and Kittel 0]. The analysis of its physical content was 
performed based on a Hubbard Hamiltonian applied to 
a dimer of magnetic sites fioj ] . The microscopic origin of 
the isotropic non-Heisenberg behavior is here reanalyzed. 
A magnetic Hamiltonian is extracted at the fourth-order 
of perturbation from a Hubbard Hamiltonian of a trimer 
of magnetic sites. In comparison to previous works, the 
derivation of the Hamiltonian is not limited to two-body 
operators. As shown here-after, the resulting effective 
magnetic Hamiltonian, which contains three-body opera- 
tors, dramatically improves the treatment of the isotropic 
deviation to Heisenberg behavior with respect to the 
Hamiltonian of Eq. [TJ 

Let us call oi, b\ and a 2 , b 2 the magnetic orbitals of the 
two magnetic sites M\ and M 2 . To simplify the deriva- 
tion of the Hamiltonian, it is assumed that orbitals a 
and b belong to different irreducible symmetry represen- 
tation of the local point group symmetry. In this way, the 
model involves only two hopping integrals t a = {a\\H \a 2 ) 
and ti, = (61 1 H | £>2 ) - This simplification has no conse- 
quences on the general conclusions drawn in this let- 
ter. The HDVV Hamiltonian model space is restricted 
to products of atomic ground states. For a system with 
two unpaired electrons per magnetic center these atomic 
states are the three m s components of the triplet (T_,To, 
T+) and the corresponding products for the dimer are: 
TqTq = ||(oi6i + bia 1 )(a 2 b 2 + b 2 a 2 )\, T + T_ = \aibia 2 b 2 \ 
and T_T + = \aib\a 2 b 2 \. The subsequent second-order 
perturbation theory derivation of the HDVV Hamilto- 
nian from a Hubbard Hamiltonian including the Hund 
on-site exchange integral Jh provides the physical con- 
tent of the exchange integral J. 
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where U a and Ub are the energies of the ionic config- 
urations Ni + Ni 3+ , i.e. the Coulombic repulsion of two 
electrons in the same orbital a or b. The biccntric direct 
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exchange integrals K a±a2 and K^b? give a small ferro- 
magnetic contribution to J and are neglected. 

The analysis of the physical content of the interactions 
J eff and \ eB in Eq. Q] can be obtained from the deriva- 
tion of a magnetic Hamiltonian at the fourth-order of 
perturbation from the Hubbard Hamiltonian on a bicen- 
tric system. Omitting terms proportional to U~ 3 that 
would bring negligible contributions and make the pre- 
sentation heavier, one gets the following expressions: 
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+(Si ■ S k ) ■ (St ■ Sj) - {Si ■ Si) ■ (S, ■ S k )](5) 

Here, Jifji k = B 2^ k i n t ne derivation that only con- 
siders the non-Hund state Sq. These interactions are non- 
zero for ± J. 



t t 

where B = jj — ^- is usually non-zero. One should how- 
ever note that in some systems the overlap between the 
atomic orbitals of the two centers are equal for symmetry 
reasons, leading to t a =tb and hence B=Q if U a ~Ub- How- 
ever, in most cases one hopping integral is dominant and 
J and B are of the same order of magnitude. The largest 
contributions to the biquadratic exchange (J^- — 
are provided by configurations involving a locally excited 
non-Hund singlet state Sq: 



So 



a\bi + b\a\ 



(4) 



The energy of the configurations built from a product 
of an atomic ground state on one magnetic site and Sq 
on the other is only 2Jh higher than the energy of the 
HDVV Hamiltonian model space functions while other 
outer space configurations are higher in energy. 

The complete analysis of the magnetic Hamiltonian 
at the fourth-order of perturbation from the Hubbard 
Hamiltonian matrix of a trimer of magnetic sites is rather 
elaborate, since the corresponding Hubbard space con- 
tains 400 determinants. This derivation is, however, 
much easier if one limits the extraction to the config- 
urations external to the model space with the largest 
contribution to the fourth-order in the bicentric system, 
namely the non-Hund state So- As will be justified later 
and confirmed by the numerical illustration of the theo- 
retical analysis, the neglect of the other configurations 
has actually no significant consequences. Unlike the 
derivation performed for the two-center system, the mag- 
netic Hamiltonian is no longer restricted to two-body op- 
erators when the number of magnetic sites exceeds two. 
For the trimer, locally excited non-Hund states gener- 
ate interactions which couple the three magnetic sites 
through a three-body operator. To illustrate the mecha- 
nism involved in this operator, two pathways are sketched 
in Fig. [T]that generate a coupling between \T$T + TJ) and 
|T + T_Xb) at the fourth-order of perturbation. The sys- 
tematic evaluation of all these interactions leads to the 
following Hamiltonian: 




FIG. 1: Two pathways of the fourth-order quasidegener- 
ate perturbation theory interaction between T+T-Tq and 
TqT+T-. The coupling goes through ionic configurations, 
then a function involving a non-Hund state and finally again 
ionic configurations. The local spin function To and So are 
respectively distinghuished by the sign + and - in upperindex 
of the parenthesis. 

To compare the two Hamiltonians (with and without 
the three-body operator) and their ability to reproduce 
the deviations to a strict Hciscnberg behavior, the spec- 
trum of both Hamiltonians will be compared to the 'ex- 
act' iV-electron spectrum of two-center and three-center 
embedded clusters of the La2Ni04 perovskite. These 
spectra are determined with ab initio quantum chemi- 
cal schemes that approximate as accurate as possible the 
cigen functions of the exact electronic Hamiltonian. The 
embedded cluster method consists in a highly correlated 
treatment of the electronic structure of a fragment of 
the material. This cluster contains a limited number of 
magnetic sites (in the present case, two or three Ni ions) 
and their nearest ligands (O sites) and is embedded in 
a set of pseudo potentials and point charges that repro- 
duce the main electrostatic effects of the crystal in the 
cluster region. The zeroth-order configuration interac- 
tion space of the electronic structure description, called 
the Complete Active Space (CAS), is the Hubbard model 
space. Hence, non-dynamical correlation effects are ex- 
actly treated. Next, a variational treatment of the Differ- 
ence Dedicated Configuration Interaction (DDCI) space 
is performed to include the screening effect, responsible 
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CAS 

CAS+DDCI2 

CAS+DDCI 

CAS(ext)+DDCI2 



J min 

7.40 
19.13 
24.35 
26.99 



Dimer 

Jmax 

7.55 
20.22 
26.46 
29.52 



2.7 
5.5 
8.3 
9.0 



7.33 
18.57 



Trimer 

J max 

7.66 
21.60 



S(%) 
4.4 
15.09 



FIG. 2: Left strongly localized orbitals. 



TABLE I: Exchange integrals (in meV) extracted from the 
ab initio dimer and trimer spectra. The active space in the 
CAS(ext)+DDCI2 calculation is extended with the 2p(a) or- 
bital of the bridging oxygen. 



for the decrease of the on-site Coulomb repulsion U , and 
consequently, the increase of the exchange integral J e " . 

The DDCI method is implemented in the CASDI code 
111 ], and is one of the most precise ab initio methods 
available for the treatment of magnetic systems. The 
method has frequently been applied to extract accurate 
effective interactions [12| in strongly correlated materials. 
The magnetic orbitals are optimized with the CASSCF 
method (Complete Active Space Self Consistent Fied). 
These orbitals are strongly localized on the Ni 2+ ions 
with small derealization tails on the neighboring O, 
see Fig. [2] Due to the larger Ni-0 distance in the z- 
direction, the Ni-3d(z 2 ) orbitals are less destabilized by 
the crystal field in comparison to the Ni-3c?(a; 2 — y 2 ) or- 
bitals. For symmetry reasons, the electron occupying the 
b\ = di(z 2 ) (respectively a\ = di(x 2 —y 2 )) orbital of atom 
1 can only delocalize into the 62 = di{z 2 ') (respectively 
a 2 = d 2 (x 2 — y 2 )) orbital of atom 2. Hence, only the t a 
and tf, hopping integrals are non zero. The electronic cir- 
culation between the left and right orbitals is clearly less 
important for the d(z 2 ) orbitals than for the d(x 2 — y 2 ) 
orbitals, for which a super exchange mechanism occurs 
involving the 2p(a) orbital of the bridging oxygen. As a 
consequence the B-terms are non-zero (see Eq. [3]) and 
one may expect that the three-body operator is active in 
this system. 

The ab initio spectrum of both the dimer and a lin- 
ear trimer have been calculated at several levels of cor- 
relation; CASSCF, CAS+DDCI2, CAS+DDCI, and ex- 
tended CAS+DDCI2 [H, [3. The deviation from the 
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Heisenberg behavior may be appreciated by the spread 
of the J-valucs which are extracted from the energy 
difference between several spin states. The minimal 
Jmin and maximal J max values are reported in Table [H 
as well as the percentage of variation in J defined as 
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In the first place, we observe 



that the calculated values of J e ^ compare very well with 
the experimental value J=30 meV [15j. The calculation 
of the trimer at the CAS+DDCI level is unfortunately 
impossible due to the large size of the configuration in- 
teraction space. Therefore, the discussion of the model 
Hamiltonians will be based on the comparison with the 
CAS+DDCI2 spectrum. Although, this spectrum in- 
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TABLE II: Effective interactions (in meV) extracted from the 
CAS+DDCI2 spectrum of the trimer. e is the percentage of 
deviation accounted for by the model. 



eludes only partially the important screening effects, it 
preserves the nature of the interactions and their physical 
content. 

The interactions J e ^ (bilinear exchange), X e ^ (bi- 
quadratic exchange), and 7 e ^ (three-body interaction) 
are optimized to minimize the difference between the 
spectrum of the model Hamiltonians and the ab initio 
spectrum. Table [XT] collects the results of the parame- 
ter optimization. One must note here that the numerical 
procedure of optimization of the interactions phenomcno- 
logically takes into account contributions of any outer 
space configurations such as the di-ionic and other non- 
Hund configurations that were neglected in the derivation 
of Eqs. [3] The deviation of each state can be appreciated 
by the difference between the model and the ab initio en- 
ergy. To get an average deviation per state we have added 
all these energy differences and divided by the spectrum 
width and the number of states. A percentage of devia- 
tion to the Heisenberg behavior (e) is then defined such 
that the HDVV Hamiltonian (case (a)) reproduces 0% 
of the deviation while the ab initio spectrum reproduces 
100 %. 

The introduction of the biquadratic first neighbor ex- 
change interaction (A^T, case (b)) only accounts for 7.5% 
of the deviation observed in the ab initio spectrum. 
Fig. [3] shows that the Hamiltonian given in Eq. [1] can- 
not remove the degeneracy between the 3 B3„ and the 
5 A g states, where this loss of degeneracy is actually the 
most remarkable difference between the spectrum of the 
HDVV Hamiltonian and the ab initio spectrum. The 
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FIG. 3: Comparison of the model spectra with the 
CAS+DDCI2 spectrum, (a) Heisenberg model with J-ff , (b) 
Heisenberg model extended with Ajf , (e) magnetic Hamilto- 
nian given in Eq. [5] including , Aff and the three-body 
operator interactions 7^ 723* , (f ) magnetic Hamiltonian given 
in Eq. [5]extended with the second neighbor interactions Jif , 
A 13 , and 7^ 7^ . 

second neighbor interaction does not significantly 
improve the model spectrum (case (c) and (d)). Table 
HI! indicates that is small and e remains quite small. 
Case (e) considers the Hamiltonian given in Eq. [5] with 
bilinear and biquadratic exchange, and the three body 
interaction parametrized by 7^ ^23 ■ Initially, only first 
neighbor interactions have been optimized. The three 
body operator introduces a major improvement in the 
treatment of the deviations from the HDVV Hamiltonian 
observed in the ab initio spectrum. The model spectrum 
reproduces 85% of the deviations and is able to lift the 
degeneracy between the 3 B 3u and the 5 A ff states. 7^f 723 
is of the same order of magnitude as A^, which in turn 
has approximately the same value as in case (b). The 
increase of e is not just an effect of the increasing num- 
ber of parameters. Notice that it is the same as in case 
(d), but the differences with the ab initio spectrum are 
significantly smaller than after the fitting without the 
three-body interaction. Finally, the values of the inter- 
actions have been refined considering second neighbour 
interactions (case (f)) in order to evaluate the consistency 
of the extraction. The number of extracted parameters is 
now equal to the number of energy differences, and hence, 
the percentage of deviation e accounted for is 100%. It 
is, however, interesting to note that the values of the 
second neighbor interactions are very small and the first 
neighbor interactions almost identical to those obtained 



in case (e), thus validating the use of the smaller set of 
three parameters for a subsequent treatment of the col- 
lective properties. 

In conclusion, the results show that a biquadratic ex- 
change is a necessary ingredient but is not sufficient to 
account for deviations to strict Heisenberg behavior in 
extended systems. In the here-considered example, the 
three-body operator happens to be responsible for the 
main contribution to the deviation. Recalling that the 
three body interaction can only be active when jf- ^ ^- , 
it is worth noting that this is generally the case and 
that this interaction is most probably of importance in 
other systems. At the fourth-order of perturbation, four- 
body operators cannot appear in the Hamiltonian of a 
tetrameric cluster. The magnetic Hamiltonian given in 
Eq. [5] should therefore be accurate for the treatment 
of extended systems and strong contributions may be 
expected of the three-body operator interactions to the 
collective properties. This will be the subject of a forth- 
coming paper. 
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